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RANS  simulation  of  the  separated  flow  over  a 
bump  with  active  control 

By  Gianluca  laccarino,  Claudio  Marongiuf, 

Pietro  Catalanof,  Marcello  Amatof 


1.  Introduction  and  motivations 

Active  Flow  Control  (AFC)  is  an  attractive  technique  to  increase  the  aerodynamic 
efficiency  of  lifting  and  control  surfaces.  The  working  principle  is  to  use  localized  suction 
or  injection  to  modify  the  characteristics  of  the  boundary  layer  flow  near  regions  of 
separation,  thus  limiting  the  associated  losses.  The  effect  of  steady  suction  or  blowing 
has  been  studied  in  great  detail  in  the  past  (see  Braslow  (1999)  for  a comprehensive 
review),  especially  with  reference  to  high-lift  configurations.  The  major  limitation  to  the 
applicability  of  conventional  AFC  devices  is  the  need  to  provide  or  discharge  a constant 
supply  of  flow  for  blowing  or  suction,  respectively,  which  requires  complicated  piping, 
additional  energy  supply,  and  causes  efficiency  losses.  In  the  last  decade,  a new  AFC 
device,  namely  the  synthetic  jet  actuator,  has  been  introduced,  which  eliminates  most  of 
these  drawbacks.  In  this  system  the  net  mass  flow  is  zero,  because  a membrane  within 
a small  cavity  produces  blowing  and  suction  alternatively.  The  performance  of  synthetic 
jets  appear  to  be  extremely  encouraging,  but  most  of  the  experimental  analysis  and 
numerical  studies  are  performed  at  low  Re3molds  number.  The  current  focus  is  on  the 
application  of  such  devices  to  turbulent  flows  especially  in  the  aeronautical  industry. 

The  fluid  mechanics  of  s)mthetic  jets  at  low  Reynolds  numbers  is  well  understood 
(Glezer  & Amitay  2002)  and  their  difference  with  respect  to  conventional  (continuous) 
jet  has  been  elucidated  experimentally  (Smith  & Swift  2003).  In  particular,  it  has  been 
shown  that  the  vortex  shedding  behind  a circular  cylinder  can  be  completely  suppressed 
in  the  laminar  regime  (Glezer  k Amitay  2002)  and  substantial  drag  reductions  can  be 
obtained  at  moderate  Reynolds  numbers  (Amitay  et  al.  1997;  Catalano  et  al.  2002). 

The  objective  of  this  paper  is  to  investigate  the  accuracy  of  Reynolds- Averaged  Navier- 
Stokes  (RANS)  techniques  in  predicting  the  effect  of  steady  and  unsteady  flow  control 
devices.  This  is  part  of  a larger  effort  in  applying  numerical  simulation  tools  to  investigate 
of  the  performance  of  synthetic  jets  in  high  Reynolds  number  turbulent  flows.  RANS 
techniques  have  been  successful  in  predicting  isolated  synthetic  jets  as  reported  by  Krai 
et  al.  (1997).  Nevertheless,  due  to  the  complex,  and  inherently  unsteady  nature  of  the 
interaction  between  the  S3mthetic  jet  and  the  external  boundary  layer  flow,  it  is  not  clear 
whether  RANS  models  can  represent  the  turbulence  statistics  correctly. 

An  extensive  computational  and  experimental  investigation  of  turbulent  flow  with 
active  control  is  ongoing  at  NASA  Langley.  A workshop,  CFDVAL2004,  will  be  held 
in  March  2004  under  joint  sponsorship  by  NASA  and  ERCOFTAC  with  the  objective 
of  providing  benchmarks  for  CFD  validation  (see  cfdval2004.larc.nasa.gov).  The 
present  work  is  directed  towards  participation  in  this  workshop. 
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2.  Numerical  techniques 

2.1.  U-ZEN 

The  CIRA  U-ZEN  code  solves  the  compressible  RANS  equations  around  complex  aero- 
nautical configurations  using  multiblock  structured  grids.  The  numerical  discretization  is 
based  on  a second-order  cell-centered  finite-volume  method  with  explicitly  added  (fourth- 
order)  artificial  dissipation.  The  unsteady  solution  procedure  is  based  on  the  dual  time 
stepping  method  where  a pseudo  steady-state  problem  is  solved  at  each  time  step.  Con- 
ventional convergence  accelerators,  including  geometrical  multigrid  and  residual  smooth- 
ing, are  used.  Several  turbulence  models  are  available  in  U-ZEN:  for  the  numerical  sim- 
ulations presented  in  this  paper  only  the  one-equation  Spalart-Allmarass  (SA)  model 
(Spalart  &:  Allmaras  1994)  and  the  two-equation  k-co  (KW)  model  (Wilcox  1993)  have 
been  used. 


2.2.  FLUENT 

FLUENT  is  a commercial  CFD  code  that  solves  the  RANS  equations  on  hybrid  unstruc- 
tured grids.  It  uses  a second-order  upwind  discretization  based  on  the  SIMPLE  pressure- 
velocity  coupling  and  the  formulation  can  accommodate  compressible  flows.  Dual  time 
stepping  is  used  to  obtain  time  accurate  simulations  and  an  algebraic  multigrid  tech- 
nique is  used  to  accelerate  convergence  within  each  time  step.  A multitude  of  turbulence 
models  and  variants  are  available  in  Fluent.  In  this  work  the  SA  model  and  Durbin’s 
V2F  four-equation  model  (Durbin  1995)  are  used.  The  latter  model  was  implemented  in 
FLUENT  using  the  User  Subroutines  (laccarino  2001). 


3.  Preliminary  validation 

As  a preliminary  step  toward  applying  the  two  codes  to  AFC  problems,  flows  around 
bluff  bodies  at  low  Reynolds  number  are  considered.  In  particular,  FLUENT  has  been 
successfully  applied  to  compute  the  vortex  shedding  in  two-  and  three-dimensional  flows 
(Ooi  et.  al.  2002);  the  same  problems  have  been  carried  out  using  U-ZEN.  Only  a few 
results  are  reported  for  the  flow  around  the  circular  cylinder  at  two  Reynolds  numbers 
{Ren  = 100  and  Rcd  — 3,900)  corresponding  to  unsteady  laminar  and  turbulent  flows. 
The  results  are  summarized  in  the  Table  I. 


U-ZEN 

FLUENT  Exp.i 

St 

0.164 

0.143 

0.165 

Cd 

1.373 

1.377 

1.340 

Cl 

0.292 

0.340 

0.325 

Lr 

1.450 

1.380 

1.460 

^min 

-0.180 

-0.185 

-0.180 

Rep  — 100 
* Zdravkovich  1997 


U-ZEN 

FLUENT 

LESi 

Exp.2 

St 

0.208 

0.233 

0.203 

0.215 

Cd 

0.865 

1.050 

1.000 

0.980 

Lr 

1.620 

1.360 

1.360 

1.330 

^min 

-0.180 

-0.185 

-0.180 

-0.180 

Rep  - 3, 900 

' Beaudan  & Moin  1994,^  Zdravkovich  1997 


Table  I.  Computational  results  for  the  flow  around  the  circular  cylinder.  St  is  the  Strouhal 
number,  Cp  and  Cp  are  the  maximum  aerodynamic  forces  during  a shedding  cycle,  Lr  the 
length  of  the  recirculation  bubble,  and  Umin  the  minimum  velocity  on  the  wake  centerline. 
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Figure  1.  Sketch  of  the  bump  studied  experimentally  by  Seifert  &c  Pack  (2002) 

U-ZEN  and  FLUENT  are  in  good  agreement  with  each  other  and  with  the  experimental 
data.  Additional  simulations  were  carried  out  for  the  flow  over  a square  cylinder,  and 
similar  agreement  was  obtained. 


4.  Turbulent  flow  over  a bump 

Three  problems  have  been  proposed  for  the  CFDVAL2004  Workshop:  the  first  two 
are  related  to  simple  synthetic  jet  flows  with  and  without  crossflow.  The  third  concerns 
the  study  of  a hump  model  with  various  steady  and  unsteady  separation  controls.  Our 
preliminary  work  has  been  focused  on  the  analysis  of  this  test  problem. 

The  geometry  with  a detailed  view  of  the  flow  control  slot  is  sketched  in  Fig.  1.  The 
hump  thickness  is  20%  and  the  chord  of  the  hump  relative  to  to  the  height  of  the  channel 
is  c/H  — 2/3  and,  therefore,  substantial  blockage  effects  are  expected. 

Experiments  were  carried  out  in  the  NASA  Langley  Transonic  Cryogenic  Tunnel 
(Seifert  & Pack  2002)  for  Reynolds  numbers  ranging  firom  2.4  to  26  millions  and  Mach 
number  of  0.25.  Due  to  discrepancies  observed  with  previous  CFD  analysis  of  this  model 
(see  Viken  et.  al.  2003)  a new  set  of  measurements  will  be  collected  and  made  available 
at  the  workshop.  In  this  work  we  will  the  original  set  of  experimental  data  for  an  initial 
assessment  of  the  predictive  capabilities  of  the  RANS  tools  employed. 

Structured  and  unstructured  computational  grids  have  been  generated.  The  unstruc- 
tured grid  (Fig.  2)  is  built  using  24  layers  of  quadrilaterals  in  the  boundary  layers  and 
paving  in  the  rest  of  the  domain.  Three  meshes  were  considered  to  evaluate  the  grid 
sensitivity  of  the  solution;  their  sizes  range  from  12,000  to  70,000  elements.  Note  that 
the  cavity  region  is  included  in  the  mesh  (using  triangular  paving)  to  study  its  effect  on 
the  solution.  The  structured  grid  is  shown  in  Fig.  3.  It  contains  16,000  cells,  and  two 
additional  coarser  versions  have  been  considered  to  evaluate  grid  dependence.  In  this  case 
the  cavity  is  not  included  and  the  jet  control  is  modeled  through  a boundary  condition 
to  be  discussed  later. 

The  first  set  of  simulations  is  aimed  at  establishing  the  accuracy  of  the  two  RANS  codes 
used  (UZEN  and  FLUENT)  for  a case  without  flow  control.  Calculations  are  carried  out 
at  Reynolds  number  of  12  million  and  Mach  0.25.  In  Fig.  4 pressure  coefficient  distribu- 
tions along  the  bump  are  reported  using  various  turbulence  models  and  by  considering 
different  heights  of  the  channel  to  study  the  effect  of  the  blockage.  The  agreement  is 
reasonable  when  the  SA  and  KW  models  are  used  in  FLUENT  and  UZEN,  respectively. 
When  the  same  turbulence  model  (SA)  is  used  the  comparison  is  not  satisfactory,  with 
UZEN  predicting  a substantially  lower  pressure.  The  cause  seems  to  be  related  to  the 
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Figure  2.  Unstructured  computational  grid  with  details  of  the  slot  region. 


Figure  3.  Structured  computational  grid.  Note  that  the  cavity  is  not  included  but  is  modeled 

as  a boundary  condition. 


incoming  boundary  layer  (note  the  peak  at  x/c  = 0 in  Fig.  4(a))  which  appears  to  be  of 
different  thickness  in  the  various  calculations.  This  was  also  acknowledged  as  a crucial 
aspect  of  the  simulations  reported  in  Viken  et.  al.  (2003). 

The  effect  of  the  blockage  plotted  in  Fig.  4(b)  shows  that  the  separation  bubble  (from 
x/c  « 0.6)  is  strongly  affected  by  the  height  of  the  channel;  however,  the  simulations  fail 
to  capture  the  level  of  the  pressure  in  the  separated  region.  This  is  consistent  with  the 
finding  in  Viken  et.  al.  (2003)  which  has  motivated  the  new  experimental  study. 

In  Fig.  5 the  streamlines  are  reported  for  the  cases  without  control  and  with  a steady 
suction  corresponding  to  a momentum  coeflScient  = —1.4%  (where  is  defined  as 
the  ratio  of  jet  momentum  to  free-stream  momentum).  The  recirculation  bubble  in  Fig. 
5(a),  estimated  using  the  skin  friction  coefficient  (not  reported),  is  A/c  = 0.59  whereas 
it  is  reduced  to  A/c  = 0.14  with  the  control  (Fig.  5(c)).  The  experimental  bubble  length 
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O Experiments  - H/c  = 1 .5 


Figure  4.  Pressure  coefficient  distribution  over  the  bump;  Re  = 16  million,  Ma  = 0.25.  Com- 
parison  between  U-ZEN  and  FLUENT  (left);  effect  of  the  height  of  the  tunnel  (right),  cfn.  Fig. 
1. 


Figure  5.  Computed  streamlines:  Re  = 16  million,  Ma  = 0.25.  No  control  (left);  no  control 
but  model  with  the  cavity  (center),  steady  suction:  — —1.4%  (right). 


is  estimated  to  be  A/c  = 0.56  for  the  case  without  control.  No  information  is  given  for 
the  case  with  control  in  Seifert  & Pack  (2002).  In  Fig.  5(b)  an  additional  simulation  is 
reported  for  the  case  without  control  but  with  the  cavity  to  evaluate  its  effect  on  the 
solution. 

One  important  aspect  of  the  present  work  is  to  establish  the  suitability  of  representing 
the  control  jets  using  boundary  conditions  instead  of  explicitly  simulating  the  flow  in  the 
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Figure  6.  Computed  streamwise  velocity  contours.  Steady  suction  modeled  as  a boundary 
condition  (top)  and  steady  suction  applied  through  the  cavity  (bottom). 


cavity.  This  is  extremely  important  for  synthetic  jets  generated  by  a deforming  mem- 
brane which  would  require  an  extremely  complicated  and  time-consuming  moving-mesh 
calculation.  In  Fig.  6 two  simulations  of  the  steady  suction  control  case  (reported  earlier) 
are  compared.  In  one  case  the  suction  is  modeled  as  a boundary  condition  whereas  in  the 
other  the  cavity  is  included  in  the  simulation.  Prom  an  experimental  point  of  view  the 
momentum  ratio  is  known  together  with  the  total  mass  flow  rate  through  the  slot: 

p,V^H 

l/2pooV2,c 

^ PocV^C 

where  the  subscripts  j and  oo  refer  to  jet  and  inflow  conditions,  respectively,  and  Mj  is 
the  measured  mass  flow  through  the  slot.  These  two  conditions  allow  the  jet  conditions 
to  be  determined: 


Poo  2Mj 
Pj  kro  AdfC 

TZ~v~ir 

The  suction  boundary  condition  is  formulated  by  assuming  a constant  density  and 
velocity  (given  by  the  expressions  above)  and  by  assuming  a zero  pressure  gradient.  The 
results  are  presented  in  Fig.  6 in  terms  of  velocity  contours  on  the  bump  and  in  the 
vicinity  of  the  suction  slot.  The  agreement  between  the  two  simulations  (with  the  cavity 
or  with  the  suction  boundary  condition)  is  very  good.  Another  option  that  has  been 
tested  for  the  specification  of  the  boundary  condition  is  to  use  a zero  gradient  condition 
for  the  velocity  corresponding  to  an  extrapolation  from  the  inside  (instead  of  a direct 
enforcement).  This  condition  will  not  guarantee  the  correct  specification  of  the  mass  flow 
and,  therefore,  the  velocity  has  to  be  properly  scaled.  The  advantage  of  this  approach 
is  that  the  velocity  is  not  constant  in  the  slot  region  but  has  a variation  that  accounts 
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Figure  7.  Pressure  coefficient  distribution  over  the  bump:  Re  = 16  million,  Ma  = 0.25.  No 
control  (left);  steady  suction;  = —1.4%  (right). 

for  the  external  flow  even  if  the  correct  suction  integral  parameters  are  specified.  The 
difference  between  the  two  boundary  conditions  described  was  found  to  be  minimal  in 
the  present  problem. 

In  Fig.  7 the  pressure  distributions  on  the  bump  with  and  without  (steady)  control 
are  reported.  Results  obtained  using  two  turbulence  models  are  compared  with  the  ex- 
perimental data.  The  V2F  model  appears  to  be  superior  to  the  SA  model  in  predicting 
the  pressure  level  in  the  recirculation  bubble  and  in  the  successive  recovery  region  (as 
expected,  see  laccarino  2001).  In  particular,  it  is  interesting  to  note  that  the  agreement 
is  somewhat  better  in  the  controlled  case. 

Finally,  in  Fig.  8,  results  for  an  unsteady  control  case  are  reported.  They  correspond 
to  an  average  momentum  coefficient  < >=  0.95%  and  a non-dimensional  frequency 

F+  = 1.6  (normalized  by  free  stream  velocity  and  half  chord).  In  this  case  experimental 
data  are  not  available  for  comparison,  but  the  present  data  axe  is  agreement  with  the  CFD 
results  reported  in  Viken  et.  al.  (2003).  It  is  interesting  to  note  that  the  envelope  (the 
maximum  and  minimum  within  a cycle)  shows  a large  variation  of  the  friction  coefficient 
in  the  region  downstream  of  the  slot  with  a very  small  area  of  recirculating  flow. 


5.  Conclusions  and  future  work 

This  report  presents  some  preliminary  results  of  steady  and  unsteady  RANS  simu- 
lations of  flows  with  active  control.  The  main  objective  is  to  establish  the  accuracy  of 
the  predictions  with  particular  emphasis  on  (a)  turbulence  modeling,  and  (b)  control 
enforcement  (boundary  condition). 

Calculations  are  carried  out  using  an  unstructured  commercial  CFD  code,  FLUENT, 
and  a multiblock  structured  aeronautical  CFD  code,  UZEN.  Several  turbulence  model 
have  been  used,  ranging  from  a one-equation  model  to  a four-equation  model.  An  explicit 
description  of  the  control  device  (a  cavity  with  a contoured  slot)  and  its  modeling  by  a 
velocity  boundary  condition  have  been  considered  and  results  have  been  compared. 
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Figure  8.  Pressure  coefficient  distribution  (left)  and  skin  firiction  (right)  over  the  bump  in  the 
unsteady  control  case;  Re  = 16  million,  Ma  = 0.25.  Solid  line  represent  time  averaged  results;  the 
grey  area  represents  the  unsteady  envelope.  These  results  are  obtained  with  the  full  configuration 
(including  the  cavity). 


The  preliminary  conclusion  is  that  turbulence  modeling  plays  a crucial  role  in  the 
accuracy  of  the  results  and,  at  least  for  a steady  suction  control,  the  control  device  can 
be  reasonably  modeled  by  a boundary  condition. 

Future  work  will  be  focused  on  the  comparison  of  the  present  calculation  (and  ad- 
ditional simulations  in  different  conditions)  to  the  new  set  of  experimental  data  being 
provided  at  the  CFDVAL2004  workshop  and  on  the  assessment  of  the  suitability  of  mod- 
eling unsteady  control  (synthetic  jet)  by  boundary  conditions. 
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